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We present a semi-analytical treatment of both the elastic and inelastic coUisional properties near 
a p-wave Feshbach resonance. Our model is based on a simple three channel system that reproduces 
more elaborate coupled-channel calculations. We stress the main diflterences between s-wave and 
p-wave scattering. We show in particular that, for elastic and inelastic scattering close to a p-wave 
Feshbach resonance, resonant processes dominate over the low-energy behaviour. 
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INTRODUCTION 

The observation of molecular gaseous Bose-Einstein 
condensates (BEG) and the subsequent experimental 
study of the BEC-BCS crossover [llill0,'5^ were made 
possible by the possibility of tuning interatomic inter- 
actions using a magnetic field (the so-called Feshbach 
resonances). Although all these experiments were based 
on s-wave interatomic interactions, it is known from con- 
densed matter physics that superfluidity of fermionic sys- 
tems can also arise through higher order partial waves. 
The most famous examples of this non-conventionnal su- 
perfluidity are ^He @, for which the Cooper pairs spawn 
from p-wave interactions, and high-Tc superconductivity, 
in which pairs are known to possess d-wave symmetry [3 ■ 
Recent interest in p-wave interactions in cold atom gases 
stemmed from these possibilities and resulted in the ob- 
servation of p-wave Feshbach resonances in and 
^Li 0, 01 , as well as theoretical studies on the super- 
fluidity of cold atom interacting through p-wave pairing 

Iii,i2.j. 

The present paper is devoted to the study of p-wave 
interactions close to a Feshbach resonance and it derives 
some results presented in 0. In a first part, we present 
the model we use to describe both elastic and inelastic 
processes that are discussed in the second part. We stress 
in particular the main qualitative differences between p- 
wave and s-wave physics and show that contrarily to the 
case of s-wave that is dominated by low energy physics, 
p-wave scattering is dominated by a resonance peak as- 
sociated to the quasi bound molecular state. Finally, 
we compare our analytical results to numerical coupled- 
channel calculations. 



MODEL FOR P-WAVE INTERACTIONS 

We consider the scattering of two identical particles of 
mass m. As usual when treating a two-body problem. 



we work in the center of mass frame and only consider 
the motion of a fictitious particle of mass m/2 interact- 
ing with a static potential. In order to study the p-wave 
Feshbach resonance, we use a model based on the separa- 
tion of open and closed channels. In this framework, the 
Feshbach resonance arises in an open channel as a result 
of the coupling with a closed channel lU]- At resonance, 
scattering properties are dominated by resonant effects 
and we can neglect all "background" scattering (ie. we 
assume there is no scattering far from resonance). 

1. We restrict ourselves to a 3-channel system, la- 
belled 1,2 and 3 which correspond to the different 
two-body spin configurations (Fig. Channels 
|1) and 1 2) are open channels. We focus on the 
situation were atoms are prepared initially in state 

Atoms may be transferred to state |2) after an 
inelastic process. Channel |3) is closed and hosts 
the bound state leading to the Feshbach resonance. 

Let us consider for instance the case of ^Li atoms 
prepared in a mixture of \F = 1/2, mp = 1/2) 
and \F' = 1/2, = -1/2)- In this system, the 
only two-body decay channel is associated with the 
flipping of a m'p = —1/2 atom to m'p = 1/2. 
If we denote by {mF,'m'p) the symmetrized lin- 
ear combination of the states \F — 1/2, mp) and 
\F' = 1/2, m'p), then |1) = (1/2,-1/2) and |2) = 
(1/2,1/2). 

2. The Feshbach resonances studied here are all lo- 
cated at values of the magnetic field where the 
Zeeman splitting is much larger than the hyperflne 
structure. In flrst approximation we can therefore 
assume that the internal states are described by 
uncoupled electronic and nuclear spin states. In 
the absence of any dipolar or hyperflne coupling 
between the electronic singlet and triplet mani- 
folds, we assume we have no direct interaction in 
channels 1 and 2 so that the eigenstates are plane 
waves characterized by their relative wave- vector k 
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and their energy Ei{k) = h'^k'^/m (channel 1) and 
£'2(fc) = -A + n^k'^lm (channel 2). A > is the 
energy released in an inelastic process leading from 
1 to 2. A can be considered as independent of the 
magnetic field and is assumed to be much larger 
than any other energy scales (in the case relevant 
to our experiments, A/h ^ 80 MHz is the hyperfine 
splitting of ®Li at high field) . 

3. In channel 3, we consider only a p-wave bound state 
nesting at an energy d quasi-resonant with channel 
1. In the case of ^Li atoms in the = 1/2 hy- 
perfine state, S = 2^b{B — Bq), where B is the 
magnetic field and Bq is the position of the "bare" 
Feshbach resonance. If the projection of the an- 
gular momentum (in unit of H) is denoted by mu 
for a quantization axis chosen along some vector u, 
the eigenf unctions associated to this bound state 
can be written as g{r)Y{^"{9,(j)), where {r,9,<l)) is 
the set of polar coordinates and the Y™ are the 
spherical harmonics. 

4. The coupling V between the various channels af- 
fects only the spin degrees of freedom. Therefore 
the orbital angular momentum is conserved during 
the scattering process and we restrict our analy- 
sis to the p-wave manifold. This is incontrast to 
the situation in heavy alkalis where incoming par- 
ticles in s-wave can be coupled to molecular states 
of higher orbital angular momentum |l4l ll^l . 

We assume also that the only non-vanishing ma- 
trix elements are between the closed and the open 
channels (ie. (1,2|T>|3) and (3|y|l,2)). 

Let us denote a state of the system by |q;,x), where 
a G {1,2,3} and x describe respectively the internal 
(spin) and the orbital degrees of freedom. According to 
assumption Q), the matrix clement (a, x|1^|q;', x') is sim- 
ply given by: 




FIG. 1: The p-wave model: we consider three internal states, 
labelled jl), |2) and |3). States |1) and |2) are two open chan- 
nels corresponding respectively to the incoming and decay 
channels. The released energy in an inelastic collision bring- 
ing atom initially in |1) to |2) is denoted A. State |3) is a 
closed channel that possesses a p-wave bound state of energy 
S nearly resonant with state Finally, we assume that these 
three channels interact through a potential V acting only on 
the internal states and coupling the two open channels to the 
closed channels. 
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(a,k|l/|3,mk) =i(5™,,o^^(a|1^|3) J g* {r)n{kr)T^ Ar, 

(2) 

where is a quantization volume. Since for small k we 
have ji{kr) ~ kr/i, the matrix element (a, k|y |3, Wk) 
takes the general form 



(Q;,k|y|3,mk) 



(3) 



.x\V\a',x')^{x\x'){a\V[ 



(1) 



where Fa(k) has a finite (in general non zero) limit when 
k goes to zero. 

Later on, we shall also need the coupling between |a, k) 
and |3,mk' = 0) (that will be denoted by |3,0k')), where 
the momentum k and the direction of quantization k' 
are no longer parallel. The calculation presented above 
yields readily 



and is therefore simply proportional to the overlap (xlx') 
between the external states. 

Let us now particularize to the case where a G 
{1,2}, and |x) = |k) is associated to a plane wave 
of relative momentum ?ik. According to hypothesis 
this state is coupled only to the closed channel 
|3). Moreover, using the well known formula e^^^ — 
Y,ii^\/^^WTT)ji{kr)Yi°{e,(i)), where the ji are the 
spherical Bessel functions, we see that |x) is coupled only 
to the state \a' = 3,k = 0) describing the pair in the 
bound state \a' = 3) with zero angular momentum in 
the k direction. The matrix element then reads 



(a,k|y|3,0k') = ^^^(OklOkO = ^^^cos(k,k'), 



where (k, k') is the angle between k and k' [l6j |. 



(4) 



T-MATRIX 

From general quantum theory, it is known that the 
scattering properties of a system are given by the so- 
called T-matrix T. It can be shown in particular that 
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\a,k> |ari,A:i> |Q'2.*2> k«.*n> \a'Jc'> 



FIG. 2: Diagrammatic expansion of the T-matrix. The full lines (resp. dashed) represent free atoms (resp. molecules). |a, k) 
is the scattering state of the two particles, in the internal state a — 1,2. |3, Ok) represents the state of a p-wave molecule with 
orbital angular momentum component zero on k direction. 



T is given by the following expansion in power of the 
coupling potential 



T(E) = V + VGq{E)V + VGo{E)VGo{E)V 



(5) 



From Eqn. (O , we see that replacing Fa (q) by its value 
at (7 = leads to a divergence. This divergence can be 
regularized by the use of counter terms in the integral, 
namely by writing that 



where GoiE) = 1/{E - Hq) emd Hq = H - V is the free 
hamiltonian of the system. 

Let us consider |a,k) and |Q;',k'), two states of 
the open channels and let us set T^a' (k, k' , -B) = 
(a, k|T(i?)|a', k'). According to formula this matrix 
element is the sum of terms that can be represented by 
the diagram of Fig. |21 and we get after a straightforward 
calculation 



Taa'{-k,-k',E) = ^Fa{k)F:^,{k')Y,RnHErGi"'\Er+' 



\FJl)\ 



mq 



E- 

dq 



- (E-E^m 



-riE-Eaio)) 

n 



dq 
6^ 



67r2 



where we have assumed that Fa was decreasing fast 
enough at large q to ensure the convergence of the in- 
tegrals. |FQ(g)p can now be safely replaced by Aq = 
|Fa(0)p in the first integral and we finally get 



ri=0 



Here Gq™^ {E) = 1/{E — 5) is the free propagator for 
the molecule, A = Ai + A2 with 



. Aq m / m 



with 



EaiO)) 



3/2 



-So,a-Vo.{E-Ea{0)), 



AaiE) 



q'dq \Faiq)\^ 

{2TlfE-Ea{q) 



(6) 



results of the integration on the loops, and finally 



^ri = / cos(k, ki) cos(ki, k2)... cos(k„, k'). 



where Vtp is the solid angle associated to kp, arises from 
the pair breaking vertices |3,0ki) |ai+i,ki+i). This 
last integral can be calculated recursively and we get 
i?„ = (47r/3)"cos(krk'); that is for the T-matrix 



T ,fk k' E) 1 kk'Fa{k)F:,{k') ^ 

with Sq, — ATiKa/i. 

This expression can be further simplified since, accord- 
ing to Eqn. (O, the width of Fa{q) is of the order of l/i?e, 
where Re is the characteristic size of the resonant bound 
state. In the low temperature limit, we can therefore 
expand Eq with the small parameter kRf,. 



\F.{q)\ 



2 mq^ dq 



If we assume that the release energy A is much larger 
ran E and if we set ^ 
get for the T-matrix 



than E and if we set do = Sq^^ + 6^^^ and f] = rji + rj2, we 



Taa'{k,k',E) 



with 



1 fcfc'TQ(0)T*,(0)cos(k,k')/(l + 77) 



L3 



E-S + ih-f{E)/2 



(7) 



my/' (A2A3/2 + Aii;3/2) 
n^J 37r(l + ??) 

{S - So) 

1+7] ' 



We note that this expression for the T-matrix is con- 
sistent with the general theory of multichannel scattering 
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resonances ; where resonantly enhanced transitions 
to other channels are readily included. In a similar con- 
text of two open channels and a Feshbach resonance, a 
recent experiment was analyzed [isj . that involved the 
decay of a molecular state formed from a Bose-Einstein 
condensate. 



S-WAVE VS P-WAVE 

This section is devoted to the discussion of the expres- 
sion found for the T-matrix. In addition to the scat- 
tering cross section, the study of the T-matrix yields 
important informations on the structure of the dressed 
molecular state underlying the Feshbach resonance and 
we will demonstrate important qualitative differences be- 
tween the behaviours of p-wave and s-wave resonances. 

Molecular state. The binding energy Et, of the 
molecule is given by the pole of T. In the limit S ^ 6q, it 
is therefore given by: 



On the other hand, we can also write Ei, = 
{■ipin{B)\H{B)\'ilj^{B)). Since in the absence of any de- 
cay channel, the molecular state is the ground state of 
the two-body system, we can write using the Hellman- 
Feynman relation 



AMoff = = -(V'm(-S)l \^m{B)). 

In our model, the only term of the hamiltonian depend- 
ing on the magnetic field is the energy S = 2^b(^ — -So) 
of the bare molecular state in the closed channel and we 
finally have 



^Mcff — — 2^B(closed|closed). 



(9) 



If we compare Eqn. © and lO, we see that the proba- 
bility Pcioscd = (closed I closed) to be in the closed channel 
is given by 



We see that the real part of Ei, (the "physical" bind- 
ing energy) is ~ (5 and therefore scales linearly with the 
detuning S — Sq. This scaling is very different from what 
happens for s-wave processes where we expect a {S — Sq)'^ 
behavior. This difference is in practice very important: 
indeed, the molecules can be trapped after their forma- 
tion only if their binding energy is smaller than the trap 
depth. The scaling we get for the p-wave molecules means 
that the binding energy increases much faster when we 
increase the detuning than what we obtain for s-wave 
molecules (this feature was already pointed out in |0). 
Hence, p-wave molecules must be looked for only in the 
close vicinity of the Feshbach resonance - for instance, 
for 77 <C 1 (relevant for ^Li, as we show below) and a 
trap depth of 100 /xK, the maximum detuning at which 
molecules can be trapped is ~ 0.1 G. 

This asymptotic behavior of the binding energy is 
closely related to the internal structure of the molecule. 
Indeed, the molecular wave function \^,^(B)) can be writ- 
ten as a sum |open) -t- jclosed) of its projections on the 
closed and open channels, that correspond respectively 
to short and long range molecular states. If we neglect 
decay processes by setting A2 = 0, we can define the mag- 
netic moment of the molecule (relative to that of the free 
atom pair) by 



AMeff(S) = - 



dEb 35 



dB dB' 
that is, in the case of ^Li where 6 = 2^b{B — Bq) 



^l^es{B) 



(8) 



inclosed = 1/(1 + 77). 

In other words, unless r] = 00, there is always a finite 
fraction of the wave-function in the tightly bound state. 
In practice, we will see that in the case of ^Li, rj ^ 1. 
This means that the molecular states that are nucleated 
close to a Feshbach resonance are essentially short range 
molecules. On the contrary, for s-wave molecules, E^ oc 
((5 — (So)^ leads to A^off = —2^b (closed] closed) oc {S — 60). 
This scaling leads to a zero probability of occupying the 
bare molecular state near a s-wave resonance. 

We can illustrate this different behaviours in the sim- 
plified picture of Fig. |2| For small detunings around 
threshold, the p-wave potential barrier provides a large 
forbidden region, that confines the bound state behind 
this barrier. The bound state wavefunction decays ex- 
ponentially inside the barrier and the tunneling remains 
nearly negligible. Since there is no significant difference 
for the shape of the p-wave bound state for small posi- 
tive and negative detunings, the linear dependence of the 
closed channel with magnetic field will be conserved for 
the bound state, and therefore the binding energy will 
also linearly approach the threshold. We note that the 
linear dependence close to threshold can also be found 
from the general Breit-Wigner expression for a resonance, 
in combination with the p-wave threshold behavior of the 
phase shift fH |. 

The imaginary part of Ej, corresponds to the lifetime 
of the molecule. In the case of s-wave molecules for which 
two-body decay is forbidden |^ , the only source of 
instability is the coupling with the continuum of the in- 
coming channel that leads to a spontaneous decay when 
^ > (see Fig. |31a). By contrast, we get a finite life- 
time in p-wave even at S < 0: due to dipolar relaxation 
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(a) 



(b) 




FIG. 3: Effect of the centrifugal barrier on tlie bound state in 
p-wave Fesfibacfi resonances, (a) Case of a s-wave scattering: 
tlie bare molecular state goes from 5 < (full line) to 5 > 
(dashed line). In the process, the molecular state becomes 
unstable and the wave function becomes unbounded, (b) In 
the case of p-wave bound state, the presence of the centrifugal 
barrier smoothes the transition from 5 > to 5 < 0. Even for 
5 > 0, the wave- function stays located close to the bottom of 
the well. 



between its constituents, the molecule can indeed spon- 
taneously decay towards state |2). For (5 0, the decay 
rate 70 associated to this process is given by: 



70 = 7(0) = 



A, 



1 + V iT^lT- 



3/2 



Elastic scattering. The scattering amplitude /(k, k') 
for atoms colliding in the channel 1 can be extracted from 
the T-matrix using the relation 



/(k,k') 



mL 
4^ 



rii(k, k',E^ h^k^/m) 



that is 



/(k,k') = 



mAi / A:2cos(k,k')/(l + ?7) 
Anh^ \h^k^/m-5-ih-f/2 



The cos(k, k') dependence is characteristic of p-wave 
processes and, once again, this expression shows dramatic 



differences with the s-wave behavior. First, at low fc, 
/(k, k') vanishes like k'^. If we introduce the so-called 
scattering volume Vg defined by 



then we have 



f{ky) = -Vsk^ 



-m\i 



-m\i 



if we neglect the spontaneous decay of the molecule. We 
see that in this approximation, the binding energy E^, of 
the molecule is given by 



mAi 



1 



In s-wave processes, the binding energy and the scat- 
tering length are related through the universal formula 
iJb = —It? /ma?. This relationship is of great impor- 
tance since it allows to describe both scattering prop- 
erties and the molecular state with the sole scattering 
length, without having to care with any other detail of 
the interatomic potential. In the case of p-wave, we see 
that no such universal relation exists between the scat- 
tering volume and i?b , a consequence of the fact that we 
have to deal with short range molecular states, even at 
resonance. We therefore need two independent parame- 
ters to describe both the bound states and the scattering 
properties. 

In the general case, the elastic cross-section (ToI is pro- 
portional to l/p. According to our calculation, we can 
put CTei under the general form 



(10) 



where E = T?k'^ /m is the kinetic energy of the relative 
motion and C is a constant depending on the microscopic 
details of the system. Noticeably, the energy dependence 
of the cross-section exhibits a resonant behavior dX E = 5 
as well as a plateau when E —f 00, two features that were 
observed in the numerical coupled channel calculations 
presented in 0| . Once again, this leads to physical pro- 
cesses very different from what is expected in s-wave scat- 
tering. Indeed, we know that in s-wave, we have / ~ —a, 
as long as A:a <C 1. Since a is in general non zero, the low 
energy behavior gives a non negligible contribution to the 
scattering processes. By contrast, we have just seen that 
in the case of p-wave, the low energy contribution was 
vanishingly small (cTei ~ E'^) so that the scattering will 
be dominated by the resonant peak E ^ S. 
Inelastic scattering. For two particules colliding 
in channel 1 with an energy E = r?k'^/m, the 
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2 3 4 5 

Energy [|jK] 



FIG. 4: Energy dependence of the elastic cross section. Dots; 
numerical closed channel calculation. The left peak corre- 
sponds to m; = and the right peak to mj = ±1. Full line: 
Fits using Eqn. IIIH 



probability to decay to channel 2 is proportional to 
/92(fc')l^i2(k, k', where p2 is the density of state in 
channel 2. Since k' is given by the energy conservation 
condition fi'k''^ jm — Is. = E, and in practice A 3> -E, 

we see that k' 



' m^jfi is therefore a constant. Us- 
ing this approximation, we can write the 2-body loss rate 
gi(E) for particles of energy E as: 



52 (i?) = 



BE 



(£;-^)2 + ;iV/4' 



where D is a constant encapsulating the microscopic de- 
tails of the potential. 



COMPARISON WITH COUPLED-CHANNEL 
CALCULATION 

The quantities such as C, £), 70 etc. that were intro- 
duced in the previous section were only phenomenological 
parameters to which we need to attribute some value to 
be able to perform any comparison with the experiment. 
These data are provided by ah initio numerical calcula- 
tions using the coupled channel scheme described in . 
The result of this calculation for the elastic scattering 
cross-section is presented in Fig. 0] The most striking 
feature of this figure is that it displays two peaks instead 
of one, as predicted by Eqn. H1U|) . This difference can 
be easily understood by noting that the dipolar inter- 
action that couples the molecular state to the outgoing 
channel provides a "spin-orbit" coupling that modifies 
the relative orbital angular momentum of the pair [l^ . 
In other word, each resonance corresponds to a different 
value of the relative angular momentum to/ (the m; = -f-l 



and TO/ = — 1 resonances are superimposed because the 
frequency shift induced by the dipolar coupling is pro- 
portional to m^, as noted in 

As the spin-orbit coupling is not included in our sim- 
plified three-level model, we take the multiple peak struc- 
ture into account by fitting the data of Fig. using a 
sum of three laws (I10|) with a different set of phenomeno- 
logical parameters for each value of the angular momen- 
tum: 



-1-1 



CmiE 



Til {E^5r,,,)^ + n'^i 



/4' 



(11) 



where is related to the magnetic field through (5m, — 
H{B - BF,7ni) and 7™, = 7o,mi + a^iE^/'^. Using this 
law, we obtain a perfect fit to the elastic as well as inelas- 
tic data obtained from the coupled channel calculations. 
The values obtained for the different phenomenological 
parameters are presented on Tabled 

From these data, we see first that the "elastic" proper- 
ties are independent of to/ . This comes from the fact that 
the elastic scattering is mainly a consequence of the hy- 
perfine coupling that does not act on the center of mass 
motion of the atoms. However, we see that both the in- 
elastic collision rate constant D and the molecule lifetime 
70 exhibit large variations with the relative angular mo- 
mentum • First, the spontaneous decay rate 70 of a 
molecule in to/ = 0,-1-1 is always larger than ~ lO^s"^, 
which corresponds to a maximum lifetime of about 10 ms. 
Second, we observe a strong reduction of the losses in the 
TO/ = — 1 channel, in which no significant spontaneous de- 
cay could be found. An estimate of 70 can nevertheless 
be obtained by noting that, since the elastic parameters 
are independent of to/, the ratio D/70 should not depend 
on TO/ (this can be checked by comparing the ratios D/jq 
for nil — and to/ = in the (-1/2,-1/2) channel). Us- 
ing this assumption we find that 70 '--^ 4 x 10"'^ s~^ both 
in (1/2,-1/2) and (-1/2,-1/2). The reason for this strong 
increase of the lifetime of the molecules in to; = —1 is 
probably due to the fact that due to angular momentum 
conservation the outgoing pair is expected to occupy a 
state with I = 3 after an inelastic process. Indeed, if we 
start in a two-body state (mF^rn'p) and if the dipolar 
relaxation flips the spin m'p, then the atom pair ends up 
in a state {mF,'m'p + 1). This increase of the total spin 
of the pair must be compensated by a decrease of the rel- 
ative angular momentum. Therefore, if the molecule was 
associated to a relative angular momentum to/ , it should 
end up with to/ — 1. In the case of to; = —1, this means 
that the final value of the relative angular momentum is 
TO/ — —2, ie. I > 2. But, according to selection rules 
associated to spin-spin coupling, the dipolar interaction 
can only change ^ by or 2. Therefore, starting from a 
p-wave {I = 1) compound, this can only lead to Z = 3. 
Let us now assume that the coupling between the molec- 
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channel 


(1/2,-1/2) 


(-1/2,-1/2) 


mi 


-1 





1 


-1 





1 


C(10-i^ cm^) 


0.22 


0.22 


0.22 


0.87 


0.88 


0.87 


D(10"" cmVK/s) 


0.00002 


0.59 


0.56 


3 X 10"^ 


1.54 


5.72 


70 (s"^) 


< 10"^ 


110 


110 


< 1 


220 


830 




0.0017 


0.0017 


0.0017 


0.0024 


0.0024 


0.0024 


V 


0.22 


0.22 


0.22 


0.23 


0.23 


0.23 


SBf (G) 


-0.0036 





-0.0036 


-0.012 





-0.012 



TABLE I: Values of the phenomenological parameters obtained after a fit to the coupled channel calculations data of Fig. |1] 
SBp is the shift between the m; — ±1 and mi — resonances. 



ular state and the outgoing channel is still proportional 
to the overlap between the two states (see Eqn. QJ, even 
in the presence of a dipolar coupling: the argument above 
indicates that the ratio p — 7o,mi=— i /lo.mi^-i between 
the decay rate of molecules in m/ = —1 and the one of 
molecules in m[ 7^ 1 is then of the order of 



P' 



J 9*{r)jz{kr)r'^dr 



J 9* {kr)r^dr 



where k — ymA/Ti"^ is the relative momentum of the 

atoms after the decay. For lithium, we have i?e ~ 3 A 
[i^j l which yields kRe ~ 7x 10~^. This permits to approx- 
imate the spherical Bessel function ji by their expansion 
at low k, ji{kr) ~ (fcr)', that is 



J g*{r)r^dr 



J g*{r)r^dr 



{kReT 



With the numerical value obtained for fci?e, we get p ^ 
2 X 10~^, which is, qualitatively, in agreement with the 
numerical coupled channels calculations presented above. 



TEMPERATURE AVERAGING 

In realistic conditions, the two body loss rate G2 needs 
to be averaged over the thermal distribution of atoms. G2 
is therefore simply given by 



G2[E) 



A{k^Tf 



/ g2{E)e-^/''^^E'/^dE. 
Jq 



The evolution of G2 vs detuning is displayed in Fig. [S] 
and shows a strongly asymmetric profile that was already 
noticed in previous theoretical and experimental papers 

This feature can readily be explained by noting that 
in situations relevant to experiments, 70 is small with 
respect to temperature. We can therefore replace 92 by 



a sum of Dirac functions centered on (5o,m, • When the 
So.rrii are positive, G2 takes the simplified form 



G2 =4V^^ 



3/2 



knT 



Moreover, if we neglect the lift of degeneracy due to 
the dipolar interaction coupling and we assume all the 
to be equal to some S , we get 



3/2 



(12) 



with D = and D/^ = E™, ^™,/7™, 0- 

For ^ < 0, and \5\ ^ k^T, we can replace the denom- 
inator of 92 by S'^ and we get the asymptotic form for 
G2 



G2 



0. ^ x: 



mi mi 



(13) 



Let us now comment the two equations (|12|l and l|13|l . 

1. We note that the maximum value of G2 is obtained 
for 6/kBT = 3/2. It means that when tuning the 
magnetic field (ie., S), the maximum losses are not 
obtained at the resonance 6 = 0, but at a higher 
field, corresponding to ^ = 3kBT/2. For a typical 
experimental temperature T = 10 /iK, this corre- 
sponds to a shift of about 0.1 G. 

2. Similarly, the width of 6*2 (^) scales like fceT"- Ex- 
pressed in term of magnetic field, this corresponds 
to - 0.1 G for T = lO^iK. This width is a conse- 
quence of the resonance nature of the scattering in 
p-wave processes. As seen earlier, both elastic and 
inelastic collisions are more favorable when the rel- 
ative energy E = 6. When (5 < 0, the resonance 
conditions cannot be fulfilled, since there are no 
state in the incoming channel with negative energy. 
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FIG. 5: Full line: numerical calculation of the loss-rate for 
T = 10 /mK. Dotted line: asymptotic expansion 11211 . 



The scattering is then formally analogous to opti- 
cal pumping or other second-order processes and 
yields the yjP obtained in lO. When I > k^T , 
the resonance condition is fulfilled by states that 
are not populated (since for a thermal distribution, 
we populate states up to £' ~ A:bT). 



CONCLUSION 

In this paper, we have developed a simple model cap- 
turing the main scattering properties close to a p-wave 
Feshbach resonance. The analytical formulas we ob- 
tained show very good agreement with both numerical 
coupled channel calculations and experimental measure- 
ments from our group and from the Innsbruck group 
[isf. We have shown that the line shape of the reso- 
nance is very different from what is expected for a s-wave 
process: while s-wave scattering is mainly dominated by 
low energy processes, p-wave scattering is rather domi- 
nated by collisions at energies equal to that of the molec- 
ular state. Regarding p-wave molecules, we have seen 
that at resonance their wave-function was dominated by 
the short range bare molecular state. Finally, the study 
of the spontaneous decay of these molecules has shown 
a very different lifetime depending on the relative an- 
gular momentum of its constituents, since molecules in 
m; = — 1 could live 10* times longer than in m; — 0,-1-1. 
This very intriguing result might prove to be a valuable 
asset for the experimental study of p-wave molecule since 
it guaranties that m; = — 1 dimers are very stable against 
two-body losses in the absence of depolarizing collisions. 
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